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The ground state phase diagram of the spin- 1/2 Heisenberg frustrated A-chain with a fer- 
romagnetic main chain is investigated. In addition to the ferromagnetic phase, various non- 
magnetic ground states are found. If the ferromagnetic coupling between apical spins and the 
main chain is strong, this model is approximated by a spin-1 bilinear-biquadratic chain and the 
spin quadrupolar phase with spin-2 gapless excitation is realized in addition to the Haldane 
and ferromagnetic phases. In the regime where the coupling between the apical spins and the 
main chain is weak, the numerical results which suggest the possibility of a series of phase 
transitions among different nonmagnetic phases are obtained. Physical pictures of these phases 
are discussed based on the numerical results. 

KEYWORDS: A-chain, Haldane phase, spin-quadrupolar phase, Uimin-Lai-Sutherland point, ferromag- 
netism 



1. Introduction 

The interplay of frustration and quantum fluctua- 
tion has been extensively studied in a variety of low- 
dimensional quantum magnets. Even in one-dimensional 
cases, various exotic quantum phenomena such as spon- 
taneous dimerization, 1 ) l/3-plateau, 2 ~ 5 ) and noncollinear 
ferrimagnetism, 6 ' 7 ) are reported. Among them, the 
Heisenberg A-chain, which consists of a single main chain 
and apical spins interacting with the main chain as shown 
in Fig. 1, has been extensively studied as one of the sim- 
plest examples of frustrated quantum spin chains. 8 ~ 14 ) 

In the present study, we investigate the spin 1/2 
Heisenberg A-chain with a ferromagnetic main chain. 
It is assumed that one of the bonds between the api- 
cal spins and the main chain is antiferromagnetic and 
the other is ferromagnetic. These interactions introduce 
frustration. In spite of its simplicity, this model exhibits a 
variety of exotic quantum phases. It should also be noted 
that this model is closely related to other important 
theoretical models, such as a spin-1/2 ferromagnetic- 
antiferromagnetic alternating chain, 15 -' a mixed spin 
chain with frustrated side chains, 16 ) a spin-1 bilinear- 
biquadratic chain, 17 ~ 23 ) a one-dimensional Kondo neck- 
lace 24 " 29 ) and a spin 1/2 ladder with ferromagnetic 
legs. 30 " 33 ) 

In this model, the ferromagnetism of the main chain is 
destroyed by the frustrated coupling to the apical spins. 
In contrast to the mechanism of quantum destruction of 
antifcrromagnetism, which has been extensively studied 
related with the high-T c superconductivity, the mecha- 
nism of quantum destruction of ferromagnetism has been 
less extensively studied. Recently, however, experiments 
on such phenomena have been reported in one- and two- 
dimensional materials 34, 35 ) with ferromagnetic nearest 
neighbour and antiferromagnetic next nearest neighbour 
couplings. Theoretical investigation has also been car- 
ried out for corresponding models. 8, 36 ~ 38 ) However, we 
find different nonmagnetic phases near the ferromagnetic 
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Fig. 1. Structure of A-chain. 



phase in the present A-chain. 

In the next section, the model Hamiltonian is pre- 
sented. Various limiting cases are discussed in §3. The 
numerical phase diagram is presented in §4 along with 
the physical description of each phase. The last section 
is devoted to summary and discussion. 

2. Hamiltonian 

We consider the spin 1/2 Heisenberg A-chain with a 
ferromagnetic main chain represented by a Hamiltonian 



N/2 



H 



= ^(KcriTi - J F cricr l+1 - K F <Ti + iTi), (1) 



;=i 



where i"T-bonds are antiferromagnetic and Jf and K-p 
bonds are ferromagnetic. The operators eri and Ti are 
spin-1/2 operators. Periodic boundary condition is as- 
sumed unless specifically mentioned. In the following, we 
use the parametrization k = K/Jp and a = K-p/K. 
This model can also be regarded as a ferromagnetic- 
antifcrromagnetic alternating chain 15 ) with ferromag- 
netic next nearest interaction. This model is also derived 
as a limiting case of the mixed spin chains with frustrated 
side chains. 16 ) 

3. Various Limiting Cases 

3.1 0<a<<l 

For a = 0, this model can be regarded as an SU(2) 
symmetric version of the Kondo necklace 25-29 ) with a 
ferromagnetic main chain and antiferromagnetic Kondo 



2 J. Phys. Soc. Jpn. 



Full Paper 



Kazuo Hida 



interaction. Therefore, it is intuitively obvious that the 
singlet pairs sit on the if -bonds in the ground state and 
the ground state is gapful even for small k. 

This intuition can also be supported by the analytic 
results for the two-leg Heisenberg ladder with ferromag- 
netic legs and antiferromagnetic rungs for which the 
gap opens for infinitesimal rung interaction. 30 ~ 32 ) In the 
present model with a = 0, one of the leg interactions is 
absent so that the gap can develop more easily with the 
infinitesimal k. Because a finite gap cannot be destroyed 
by an infinitesimal perturbation, our model remains gap- 
ful for small a. 

3.2 < k « 1 

In this limit, the sum of the cr-spins er tot = J^^ 2 <t/ 
forms a fully polarized state with |er tot | = N/4. This 
can be regarded as a macroscopic spin with length N/4. 
Limiting the Hilbert space of cr-spins in this subspace, 
the Hamiltonian (1) reduces to 



N/2 



H, 



eff 



J-f X! ( Kcr totTl - K F <T tot Tl) 



2J F fc(l - a) 



N 



"f totTtot 



(2) 



where Ttot = Ylf^i T i an d constant terms are omitted. 
This Hamiltonian commutes with T 2 ot . Therefore, the 
ground state has |<r to t + "?"tot| = for < a < 1 (non- 
magnetic ground state) and |otot+ftot| = N/2 for a > 1 
(ferromagnetic ground state). 

Within the effective Hamiltonian, the lowest excitation 
of the effective Hamitonian (2) for < a < 1 is given by 
the state with |<r tot + r to t| = 1 and r to t = N/4, which 
has the excitation energy 



AE 



2J F fc(l - a) 
N ' 



(3) 



This vanishes in the limit N — > oo. At first sight, this 
appears to give an upper bound for the lowest excita- 
tion energy and the present system appears to be gap- 
less in the thermodynamic limit for sufficiently small k. 
However, this is in contradiction with the argument for 
small a in the preceding subsection. This contradiction 
is resolved if the gap opens in a nonperturbative way 
as exp(— const. /Vk) for small k as in the case of ferro- 
magnetic leg ladders, 30, 33 -* because the above argument 
leading to (3) is based on the perturbative argument with 
respect to k. 

With the increase in a, the ground state remains non- 
magnetic until the transition to the ferromagnetic phase 
is reached. The behavior of the present model in this 
regime is, however, nontrivial. The details will be dis- 
cussed using numerical methods in §4. 

3.3 a»l 

In this case, the spins er i+1 and n form a spin-1 fer- 
romagnetic dimer. 15 ) Within the lowest order approxi- 
mation in Jp and K, the effective Hamiltonian for these 



spin-1 operators Si = er i+1 + Ti is given by 



H. 



(i) 



N/2 



eff 



(4) 



l=i 



In the weakly frustrated regime k > 1, the effective cou- 
pling between these ferromagnetic dimers is antiferro- 
magnetic and the ground state is the Haldane phase in 
terms of Si's. On the other hand, for < k < 1, the 
effective coupling is ferromagnetic; thus the ground state 
is ferromagnetic in which all spins are aligned. Therefore, 
we may expect the Haldanc-ferromagnetic phase transi- 
tion at k = 1. However, in the close neighbourhood of the 
phase boundary k = 1 , the higher order terms come into 
play. As a result, not only the bilinear term but also the 
biquadratic term appears in the effective Hamiltonian, 
and both become of the same order of magnitude in the 
vicinity of the point k — 1 as 



H. 



(2) 



off 



N/2 ( 2\ 

J2 (jjgS l S l+1 + jjg(s l Si +1 ) (5) 



with 



J. 



(2) 
eff 



4 V a 



t(4) _ -x- 

cff " 8a' 



Jv 



To simplify the Hamiltonian, we set k = 1 in the terms of 
0{a^ 1 ) because these terms are unimportant away from 
the point k = 1. 

The exact Bethe Anszatz solution of this Hamilto- 
nian (5) is obtained for = by Uimin, 17 ) Lai, 18 ) 
and Sutherland. 19 ' In the following, this point is called 
the ULS point. At this point, the lowest spin-2 quintu- 
plet excitation degenerates with the lowest spin-1 triplet 
excitation, and both become gapless. For < J^g < 

Jgff , the spin-2 excitation is the lowest gapless exci- 
tation, 2 °-22) reflecting the spin-quadrupolar quasi-long 
range order. 23 ' This phase is called the spin-quadrupolar 

(2) 

phase. For < 0, the ground state becomes ferro- 
magnetic. Therefore, in terms of the parameters of the 
original Hamiltonian (1), the ULS point is given by 



1 



1 

2a 



(6) 



and the transition to the ferromagnetic phase takes place 
at 

1 

a' 



k = 1 



(7) 



3.4 Stability of Ferromagnetic Phase 

The single-magnon instability of the ferromagnetic 
phase can be analytically investigated. The Hilbert space 
of the state with a single magnon is spanned by the states 
\l,n), where I is the position of the flipped spin and fi = a 
or t distingushes the two types of spins. 

The Hamitonian operates on the Fourier transformed 
bases 

= "tJftsSI*'^ 6 *^* 3 *)' ( 8 ) 
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as 



AH\q,a) 



fc(l 



- cos(q) + 1 



fc(l — a exp(iq)) . . 

+ ~ o |g,T) 



Q: 



AF| ffl r> = fc(l-°exp(-w) |g>g) 



fc(l - a) 



3,t) 



(9) 



(10) 



where Ai? is the Hamiltonian measured from the energy 
of the ferromagnetic state. The eigenvalue equation reads 

AE{q) 2 - A£%)((1 - cosg) - fc(l - a)) 

-|((l-a) + fca)(l-cos«)=0 (11) 

and the stability condition E(q) > of the ferromagnetic 
phase reduces to 



k(a — 1) < 1 — cos q and a(l — fc) > 1 



(12) 



for the arbitrary q. This implies that the ferromagnetic 
ground state is stable against single-magnon excitation 
for 



a > ap = 



1-k 



(13) 



for < k < 1. The instability point ap coincides with 
the critical point (7) calculated in the limit a >> 1. For 
fc > 1, the condition (12) is never satisfied. 

The energy of the single magnon with the wave number 
q in the neighbourhood of the instability point is given 
by 



a „/ n k .„ .„ , . . 1 - COS (J 

A£ 9 = -- 1 - (1 - fc a ^p-, 

2 i_ cosg+ _L_ 



(14) 



which is negative for the arbitrary nonvanishing q for 
a < ap. Therefore, all modes become unstable at ap. 
Among them, the lowest energy mode has q = w, 
which corresponds to the short-range order in the Hal- 
danc phase. However, with the decrease in fc, the q- 
dependence of AE(q) becomes weak and AE(q) tends to 
a ^-independent value — f (1 — (1 — k)a). This implies that 
almost all modes become equally unstable at a = ap for 
small fc. This might be the origin of the complex feature 
of the phase diagram in the small fc regime presented in 
the next section. 

4. Numerical Ground State Phase Diagram 

4-1 Ferromagnetic-nonmagnetic phase boundary 

The overall ground state phase diagram is presented 
in Fig. 2. The boundary of the ferromagnetic phase is de- 
termined from the ground state level cross for the finite 
chains. This boundary has no system size dependence 
for N — 12, 16,20 and 24, and coincides with ap deter- 
mined in the preceding section. Therefore, the transition 
of the ferromagnetic phase to other phases takes place as 
a single magnon instability rather than a multimagnon 
instability. In the neighbourhood of this phase bound- 
ary, various intermediate phases are found. The details 
are discussed for large a regime and small fc regime sep- 



0.5 



Haldane phase 




Fig. 2. Overall phase diagram of present model. The boundary 
of the ferromagnetic phase is denoted by a thick curve. The de- 
tailed structure of the phase diagram drawn by thin curves are 
explained in the text. 



arately in the following subsections. 

4-2 Large a regime 

As expected from the perturbational calculation in 
§3.3, there exists a gapless spin-quadrupolar phase be- 
tween the Haldane and ferromagnetic phases for large a. 
This phase is identified as the nonmagnetic phase with 
S = 2 gapless excitations, while the Haldane phase is 
characterized by the gapped excitations with S = 1. 
The phase boundary is determined from the level cross 
of the S = 1 and 5" = 2 excitations for N = 12, 18 
and 24 chains. The system sizes are chosen by consider- 
ing the quasi-trimerized nature of the spin-quadrupolar 
phase. 20,21 ' The phase transition line extrapolated to 
N — > oo is presented in Fig. 3. Although the spin- 
quadrupolar phase obtained by our numerical method 
is quite narrow, we conclude that this phase does exist 
between the spin gap phase and the ferromagnetic phase 
because its existence is confirmed by the perturbational 
calculation in §3.3. 

To confirm the universality class of the transition be- 
tween the Haldane phase and spin-quadrupolar phase, we 
estimate the central charge c determined by the formula 

where v 8 is the spin wave velocity defined by 

v s = lim ^[E qi (N) - E g (N)]. (16) 

Here, Eg(N) is the ground state energy, E qi (N), the ex- 
citation energy of the excited state with the wave num- 
ber qi = j^j2ji an d £oo is the ground state energy per 
unit cell in the thermodynamic limit. The extrapolation 
is carried out from N = 12,18 and 24. The results are 
presented in Fig. 4, which confirms c = 2 as expected 
on the ULS point. 22 ' We have also checked c = 2 within 
the spin-quadrupolar phase as shown in Fig. 5 at a = 9. 
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Haldane phase 




Ferromagnetic phase 
Ferrimagnetic phase 



10 



a 



12 



Fig. 3. Phase diagram in large a regime. Filled squares indicate 
phase boundaries determined by the extrapolation from N = 
12, 18 and 24. Open circles and squares indicate the data for 
N = 12 and N = 24, respectively. The thick line indicates the 
single magnon instability point ctp. Thin lines arc guides for eye. 



This confirms the critical nature of this phase. It is also 
numerically checked that the lowest excited state has the 
wave number q = 2tt/3 and total spin 5* = 2 within the 
spin-quadrupolar phase, as shown in Fig. 6. 

Precise calculation suggests the presence of an ex- 
tremely narrow ferrimagnetic region with intermediate 
spontaneous magnetization between the ferromagnetic 
phase and the nonmagnetic phase. In Fig. 3, open circles 
are the phase boundaries estimated for the N — 12 chains 
and the open squares are those estimated for N = 24 
chains. The system size dependence of the ferrimagnetic- 
nonmagnetic phase boundary is also very weak. There- 
fore, we expect that this ferrimagnetic phase will survive 
in the thermodynamic limit. 

This ferrimagnetic phase can be regarded as a 
frustration-induced ferrimagnetic phase, which has been 
found in various frustrated spin chains recently 6, 7 ) This 
phase vanishes for sufficiently large a where this model 
reduces to the S = 1 bilinear-biquadratic chain, which 
has no frustration. This phase also vanishes for small k, 
where the direct ferromagnetic-nonmagnetic transition is 
expected from the argument of §3.2. Within the numeri- 
cal data, however, it is impossible to determine whether 
this phase vanishes at a finite value of a or k. 

4-3 Small k regime 

The present model has two spin- 1/2 degrees of free- 
dom in a unit cell. Therefore, it is natural to expect 
that this model undergoes Gaussian phase transitions 
between different spin gap phases as the system param- 
eters are varied. To determine the phase transition lines, 
we have carried out the numerical diagonalization with 
twisted boundary condition and evaluated the spin in- 
version parity assuming the Gaussian phase transitions 
described by the conformal field theory with the cen- 
tral charge c = 1. We twist the spins er 2 and n by an 
angle ir relative to <7i around the z axis. As discussed 
by Kitazawa 39 ) and Kitazawa and Nomura, 40 ^ the spin 




oc 12 



Fig. 4. Central charge on spin-quadrupolar-Haldane phase 
boundary. 




.89 0.9 0.91 

Fig. 5. Central charge within spin-quadrupolar phase for a = 9.0. 
0.01 



AE(q)\ a=9k=0.9 



W=12 18 24 
S=1 o □ o 
S=2 . ■ • 



« 



0.5 



Fig. 6. Dispersion relation in spin-quadrupolar phase for a = 9 
and k = 0.9 with N = 12, 18 and 24. 



inversion parity should change sign if the valence bond 
structure of the nonmagnetic phase changes. 

The difference AE tV! between the ground state energy 
with the positive parity E + and that with the negative 
parity E_ is shown in Fig. 7 for k = 0.1 and 0.2. Surpris- 
ingly, the ground state undergoes a series of phase tran- 
sitions among the states with different parities, although 
the energy difference between the two ground states is 
small for small k. The phase boundaries extrapolated to 
N — > oo are shown in Fig. 8. Let us denote the l-th tran- 
sition point by ctc\k) so that ctc\k) < ai l+1 \k). For 
< a < o^\k), the valence bonds sit on the if-bonds 
in the ground state. In this phase, the spin inversion par- 
ity with twisted boundary condition is negative because a 
valence bond sits on the twisted bond. There is no phase 
boundary between this phase and the Haldane phase in 
the large a limit. Therefore, this phase is also called the 



J. Phys. Soc. Jpn. 



Full Paper 



Kazuo Hida 5 



0.001 

AE 



tw 



• # (a) k=0A 


i i i i i 


k 




A/=24 □ 






A/=20 ■ 


0.4- 




A/=16 o 






A/=12 • 

i i i i i i i i r r f i i 









0.5 



a 1 



0.02 

A£ h „ 



0.01 



(b) /<=0.2 



A/=24 □ 
A/=20 ■ 
A/=16 o 
A/=12 • 




Fig. 7. Energy difference Ai?tw = E+ — E- of ground states of 
different parities with twisted boundary condition for (a) k = 0.1 
and (b) 0.2. 



Haldane phase in the following. 

To confirm this phase transition line, we have also 
carried out the phenomenological renormalization group 
analysis 41 -* of the scaled singlet-triplet gap NAE op with 
open boundary condition using exact diagonalization 
(N = 12, 16, 20 and 24) and density matrix renormal- 
ization group (N = 36,48 and 52). The a-depcndcncc 
of NAE op is presented in Fig. 9 for k = 0.2 and 
12 < N < 48. It should be noted that iVA£; op values 
for different system sizes cross with each other around 
the critical point a^\k = 0.2) ~ 0.74, suggesting that 
the end spins appear for a > a^K Actually, the tran- 
sition points a op calculated from the crossing point of 
NAE op coincide well with for relatively large k val- 
ues (fc>0.22) as shown in Fig 8 by small filled circles. 
For smaller values of k, however, they deviate from each 
other. The crossing points of NAE op values correspond- 
ing to Qc with I > 2 are not found. These discrepancies 
are presumably due to the logarithmic correction to E op 
resulting from the SU(2) symmetry. Therefore, we rely 
on the result of the twisted boundary method for small 
k. Neverhteless, the coincidence for /c>0.22 suggests that 
this transition is a usual Gaussian transition between dif- 
ferent valence bond ground states. From continuity, it is 
natural to expect that this will be valid even for smaller 
values of k. 

(2) 

The size dependence of ai estimated by the twisted 
boundary condition method is stronger than that of . 



(3) 

length N < 20 and is found only for N = 28. These 
features of the spin-gap phases near the ferromagnetic 
phase imply that the resonation among long ranged va- 
lence bond states 16 ^ would be required to describe these 
phases rather than a simple valence bond picture. This in 



0.2 



i i i 



Haldane phase 




Fig. 8. Phase diagram in small k regime. Large open circles in- 
dicate Gaussian phase boundaries determined by the twisted 

boundary method extrapolated from N = 16, 20 and 24 for a' 1 ^ 

(2) (3) 
and ai ' , and from N = 20, 24 and 28 for a K c > . Small open circles 

indicate otf 1 for N = 28. Small filled circles indicate the phase 
boundary determined by the phenomenological renormalization 
group analysis of open chains with 12 < N < 52, where the val- 
ues of N are multiples of 4. The thick solid line is ctp, which coin- 
cides with the numerically obtained ferromagnetic-nonmagnetic 
phase boundary. Thin lines are guides for eye. For very small k, 
numerical data are unavailabe owing to the poor convergence of 
Lanczos diagonalization. 
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Fig. 9. a-dcpcndcnce of scaled gap with open boundary condition 
for k = 0.2. 
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Fig. 10. Dispersion relations on critical points a' 1 ' and for 
k = 0.1 and k = 0.24. The system sizes are N = 12,16,20,24 
and 28. 
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2/(N 1+ N 2 ) 
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Fig. 11. System size dependence of central charge estimated for 
finite size system for (a) a = with 0.08 < k < 0.32 and 

(b)a = ai 2) with 0.1 < k < 0.2 plotted against 2/(JVi + JV 2 ), 
where Ni and JV2 are system sizes used for determining finite 
size correction to ground state energy and spin wave velocity. 
The extrapolation to (iVi + N2)/2 — ¥ 00 is carried out by fitting 
the size dependence by the third-order polynominals. 



NAE r 




Fig. 12. o-dependence of scaled gap AE P with periodic boundary 
condition for k = 0.2. 



turn suggests the possibility that further series of transi- 
tions may take place before the ferromagnetic transition 
point is reached in the thermodynamic limit. 

To confirm that these critical lines are actually de- 
scribed by the c = 1 conformal field theory, we have 
also carried out the calculation of the conformal charge 
c. Unfortunately, due to the strongly nonlinear nature 
of the dispersion relation for small wave numbers shown 
in Fig 10, the extrapolated values of c depend strongly 
on the way of extrapolation. Here, we employ the fol- 
lowing method: The coefficients of eq. (15) are deter- 
mined by fitting the ground state energies of N\ and 
A^(= N\ + 4) using eq. (15) to obtain the finite size 
estimation of A(N\, N 2 ) = v s c. The numerical data for 
the excitation energy for the same set of system sizes are 
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0.04 

0.024 
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fe=0.1 
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Fig. 13. Dispersion relation for a = 1 with k = 0.1 and k = 0.24. 
The system sizes are N = 12, 16, 20, 24 and 28. 



used to determine the finite size estimation of v s (N\, N 2 ) 
by fitting the smallest wave number excitation energy 
as £[E qi (N) - E g (N)} = v s + B/N with the constant 
B. The finite size estimation of c{N\, N 2 ) is determined 
by c(-/Vi,Ar 2 ) = A(N u N 2 )/v s (N u N 2 ). This is extrapo- 
lated to Nx,N 2 -> 00 as a function of 2/(Ni + N 2 ) as 
shown in Fig. 11. The values of c thus extrapolated seem 
to converge around unity. This result is consistent with 
the assumption that the present transition lines are de- 
scribed by the c = 1 conformal field theory, although 
other possibilities cannot be totally ruled out consider- 
ing the strong dependence on the system size and on 
the way of extrapolation. This result justifies the use of 
the twisted boundary condition method to determine the 
phase boundary. 

5. Summary and Discussion 

The ground state phase diagram of the spin- 1/2 
Heisenberg A-chain with a ferromagnetic main chain 
is determined both perturbationally and numerically. 
The overall phase diagram is divided into nonmag- 
netic, ferromagnetic and extremely narrow ferrimagnetic 
phases. For large a, the spin-quadrupolar phase with 
5 = 2 gapless excitation is shown to appear as antici- 
pated from the perturbational mapping onto the 5 = 1 
bilinear-biquadratic chain. It is remarkable that the spin- 
quadrupolar phase is found in the 5 = 1/2 model with 
only bilinear exchange interaction. This increases the 
possibility of observing this phase in real materials, al- 
though the materials described by the present model is 
not realized so far. 

In the small k regime, a series of ground states with 
different parities under twisted boundary condition are 
found near the nonmagnetic-ferromagnetic phase bound- 
ary. It is assumed that these phases are stabilized by the 
resonation of various long-distance valence bond config- 
urations. However, the ground state of this parameter 
regime is not fully understood within the present numer- 
ical studies. Here we present some arguments towards 
deeper understanding of phases in this regime. 

The energy difference A.E tw between the ground states 

with different parities is extremely small for cup > a > 

(2) 

a c (k) with small k, as shown in Fig. 7. The phenomeno- 
logical renormalization group analysis of Ai? op fails to 
find etch's with I > 2, as shown in Fig. 9. The gap with 
the periodic boundary condition is almost critical in this 
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regime as shown in Fig. 12. From these behaviors, one 
might be tempted to conclude that the ground state in 
this regime is actually gapless and Gaussian lines ob- 
tained here are stable fixed lines, that do not correspond 
to phase transitions. However, such an extended gap- 
less phase with triplet lowest excitation is unlikely in 
the present model because the gap formation term due 
to period-2 spatial modulation is always relevant at the 
SU(2) invariant point of c = 1 conformal field theory. 
The spin quadrupolar phase is also excluded because it 
should have spin-2 lowest excitation. In addition, no trace 
of quasi-trimerization is found in the dispersion relation 
within this phase as shown in Fig. 13. 

Therefore, we may speculate that the phases for ap > 

(2) 

a > a c (k) are true gapped phases with an extremely 
small gap. A similar phase is found in the spin-1/2 zigzag 
chain with ferromagnetic nearest neighbour and antifer- 
romagnetic next nearest neighbour couplings, 37 ) which 
has been fairly well investigated as a one-dimensional 
example of nonmagnetic phase that emerge from fer- 
romagnetism. 8 ' 36 ' 37 ) However, the nonmagnetic ground 
state of this model has an incommensurate spin correla- 
tion and ground state degeneracy in the thermodynamic 
limit accompanied by spontaneous dimerization. 36 ' 37 ) In 
contrast, in the present A-chain model with small k, no 
signal of ground state degeneracy is found and the dis- 
persion relation has no structure at any incommensurate 
wave numbers, as shown in Figs. 10 and 13. Presumably, 
this is because the antiferromagnetic bonds in the present 
model do not compete with each other in forming valence 
bonds. Thus, we speculate that the nonmagnetic ground 
state of the present model with small k and «<of can 
be regarded as a new nonmagnetic state with a highly 
reduced gap realized by the influence of frustrated cou- 
pling to the apical spins. Unfortunately, we do not find 
an appropriate low-energy effective field theory with a 
markedly reduced mass that describes this phase. There- 
fore, full understanding of this ground state is left for 
future studies. 

It would be an interesting experimental challenge to 
synthesize the material corresponding to the present 
model and observe the exotic nonmagnetic phases de- 
scribed in this paper. Considering that the materials with 
a A-chain geometry are already synthesized, 12 ~ 14 ) it is 
promising to expect the experimental realization of the 
present model in the near future. 
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ter, University of Tokyo. This work is partly supported 
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